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Abstract We present the spectroscopic observations for 1 1 confirmed globular clusters 
of M31 with the OMR spectrograph on 2.16m telescope at Xinglong site of National 
Astronomical Observatories, Chinese Academy of Sciences. Nine of our sample clusters 
are located in the halo of M31 and the most remote one is out to a projected radius of 
78.75 kpc from the galactic center. For all our sample clusters, we measured the Lick 
absorption-line indices and the radial velocities. It is noted that most GCs of our sample 
are distinct from the HI rotation curve of M3 1 galaxy, especially for B5 14, MCGC5, H12 
and B517, suggesting that most of our sample clusters do not have kinematic association 
with the star forming young disk of th e galaxy. We fitted th e absorption line indices with 
the updated stellar population model iThomas et al.l d2010h with two different tracks of 
Cassisi and Padova, separately, by applying the x 1 — minimization method. The fitting 
results show that all our sample clusters are older than 10 Gyr, and metal-poor (—0.91 < 
[Fe/H] < —2.38 dex). After merging the spectroscopic metallicity of our work with the 
previously published ones, we extended the cluster sample out to a projected radius of 
1 17 kpc from the galaxy's center. We found the metallicity gradient for all the confirmed 
clusters exists with a slope of —0.028 ± 0.001 dex kpc -1 . However, the slope turns to be 
—0.018 ± 0.001 dex kpc -1 for all the halo clusters, which is much shallower. If we only 
consider the outer halo clusters with r p > 25 kpc, the slope becomes —0.010 ± 0.002 
dex kpc -1 and if one cluster G001 is excluded from the outer halo sample, the slope is 
—0.004 ± 0.002 dex kpc -1 . Thus we conclude that metallicity gradient for M31 outer 
halo clusters is not significant, which agrees well with the previous findings. 

Key words: galaxies: individual (M31) — galaxies: star clusters — globular clusters: 
general — star clusters: general 



1 INTRODUCTION 

Galactic formation and evoluti onary scenarios rema in among the most important outstanding problems 
in contemporary astrophysics dPerrett et all 120021) . One way to better understand these questions is 
through detailed studies of globular clusters (GCs). These objects are often considered fossils of galac- 
tic formati on and evolution pro cesses, since they formed at the very early stages of their host galaxies' 
lifecycles dBarmbv etaUl2000t) . GCs are usually densely packed, gravitationally bound spherical sys- 
tems containing several thousands to approximately one million stars. Thus, they can be detected from 
great distances and are suitable as probes for studying the properties of extragalactic systems. 
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Located at a distanc e of approximately 780 kpc dStanek & GarnavichL Il998t iMacriL l200lb 
iMcConnachie et ail 120051) . M3 1 is the nearest and largest spiral galaxy in our Local Group. It con- 
tains a la rge number of GCs and is considered an ideal laboratory for studies of star cluste rs in external 
galaxi es. IBarmbv & Huchral d200ll) estimated the total number of GCs at 460 ± 70, while IPerina et all 
(|2010|) arrived at ^530. Both of these estimates yield much larger numbers t han for the GCs in our 
Galaxy. However, from the observational evidence collected to date (see, e.g., Rich et al~l 2005b. the 
M31 GCs and their Galactic counterparts r eveal some striking similarities dFusi Pecci et al.l 1 19941: 
iDjorgovski etafl, Il997l IBarmbv et al.1 120021) . Based on survey data from the Can ada-France-H awaii 
Telescope (CFHT) and the Wide Field Camera on the Isaac Newton Telescope (INT), lHuxorl(l2007l) con- 
cluded that M31 and the Milky Way are more similar than previously thought. However. lHammer et al.l 
(120071) compared our Galaxy and M3 1 to the local disk galaxies within the same mass range and found 
that Milky Way is an exceptional disk galaxy which did not undergo any significant merging for last 10 
Gyrs so that it lacks stellar mass, angular momentum, disk radius ad metallicity of stars in the outskirts 
while M31 is a typical disk galaxy which is shaped by relatively recent merging. This may explain why 
there are more GCs in M 3 1 than that in o ur Galaxy by a factor of 3 as the merging could lead to the 
formations of GCs. Later, lYin et al.1 d2009b found that the two galaxies are similar in the radial profiles 
of star formation rate, gas profiles and stellar metallicity distributions along the disk by studying the 
chemical evolution history of the two galaxies. The authors conclu ded that the star formation efficiency 
of M3 1 disk is twice high as that in our Galaxy. iHou et all d2009l) also compared the two disk galaxies 
and concluded that the Milky Way disk contains more gas and higher star formation rate than that of 
M31. The authors also find that the scaled abundance gradients are similar for the two galaxies. These 
recent works could provide useful clues which can explain the similarities and differences of two GC 
system. Therefore, studying the properties of the GCs in M3 1 not only improves our understanding of 
the formation and structure of our nearest large neighbor, but also of our own Galaxy. 



A large number of halo GCs in M3 1 have recentl y been discovered. T hese are important to study the 
formation history of M3 1 and its dark matter content. Huxor et ail d2004l) discovered ni ne previously un- 
known GCs in the halo of M3 1 using the INT survey. Subsequently, Huxor et ail d2005l) found three new, 
extended GCs in the halo of M31, which have characteristics between typical GCs and dwarf galaxies. 
iMackev et al.l d2006l) reported four extended, low-surface-brightness clusters in the halo of M3 1 based 
on Hubble Space TeJescope/Advanced Camera for Surveys (ACS) imaging. They are structurally very 
different from typical M31 G Cs. However, their old and metal-poor characteristics are similar to those 
of typical GCs. iHuxorl d2007l) discovered 40 new GCs in the halo of M3 1 (out to 100 kpc from the galac- 
tic center) based on INT and CFHT imaging. Some of them are also very extended. These extended star 
clusters in the M3 1 halo are very similar t o the diffuse star c lusters (DSCs) associated with early-type 
galaxies in the Virgo Cluster reported by iPeng et aTl d2006l) based on the ACS V i rgo Cl uster survey. 
However, it seems that DSCs are usually fainter than typical GCs. IMackev et al.l d2007l) reported 10 
outer-halo GCs in M31, at ~45 kpc to 100 kpc from the galactic center. Eight of these were newly dis- 
covered based on deep ACS imaging. The halo GCs in their sample are ve ry bright, compa ct, and metal 
poor, and therefore quite different from their counterparts in our Galaxy. iMa et al.l d2010h constrained 
the age, metallicity, reddening and distance modulus of B379, which is located in the halo of M3 1 , with 
the SSP model and photometry. 



In this paper, we will present our new observations on a sample of new GCs, most of them are 
located far from the galaxy center. This allows us to be able to study the properties of the M3 1 outer 
halo in more detail. The paper is organized as follows. In $2]we describe how we selected our sample of 
M31 GCs and their spatial distribution. In |J3] we reported the spectroscopic observations with 2.16 m 
telescope at Xinglong site and the data reductions from which the radial velocities and Lick line indices 
were measured. Subsequently, in Sj4] we derive the ages and metallicities of GCs with \ 2 —minimization 
fitting. We also discuss our final results on the metallicity distribution in the M31 halo. We give our 
summary in f|5] 
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2 SAMPLE SELECTION 

We selected the sources from the Revised Bologna Ca talogue of M3 1 globular c l usters and ca ndidates 
(RBC v.4, available from http://www.bo.astro.it/M31i Idalleti et all 12004 I2006L 120071 l2009h , which 
is the latest and most comprehensive M31 GC catalogue so far. It contains 2045 objects, including 
663 confirmed star clusters, 604 cluster candidates, and 778 other objects that were initially thought 
to be GCs but lat er proved to be stars , asterisms, galaxies, or Hll regions. In fact, many of the halo 
clusters were from lMackev et all (120071) . who reported 10 GCs in the outer halo of M3 1 from their deep 
ACS images, of which eight were detected for the first time (see for details in £[TJ. In our work, our 
sample clusters are completely selected from RBC v.4. We selected the confirmed and bright (< 17 
mag in V band) clusters as well as being located as far as they could from the galaxy center, where the 
local background is too bright to observe. Finally, we have 1 1 bright confirmed clusters in our sample 
and most of them are located in the halo of the galaxy. Although some of our sample clusters have 
the previous spectroscopic observations by some authors, actually those clusters lack comprehensive 
spectroscopic informations. In other words, they only have radial velocities or [Fe/H] or [a/Fe] or age 
informations. Therefore it is necessary to observe the spectra of our sample clusters systematically and 
study the ages and metallicity in detail. 

The informations of our sample GCs are listed in Table Q] including coordinates, projected radii, 
V^-band magnitudes and age estimates. All the coordinates (Cols. 2 and 3) and l/-band magnitudes 
(Col. 5) are from RB C v.4 except the V mag of EXT8 which was derived from ugriz photometry of 
IPeacock et al.l (|2010|) with the transformation equation of Ijester et al.l d2005l) as V mag of EXT8 is not 
provided in RBC v.4. The projected radii from the galaxy center r„ (C ol. 4) were calculated with M3 1 
cent er coordinate 00 : 42: 44 3 1, +41 : 16 : 09.4 dPerrett et all 120021) . PA = 38° and d istance = 785 
kpc dMcConnachie et aill2005l) . The ages (Col. 6) are from a number of previous w ork: IPeacock et al.l 
d20ld) by using the SPSS and 2MASS photometric color s: ICaldwell et al.l d2009l) by using the 6.5m 



MMT Hectospec spetra line indices and HST CMD fittings. iGalleti et aL ( 20051) bv comparing the lines 
indices with the prediction models. 



Table 1 The parameters of our sample GCs. 



ID 


R.A. 


Dec. 




V 


age 


references for ages 6 




(J2000) 


(J2000) 


(kpc) 


(mag) 






MCGC2 


00:29:44.90 


+41:13:09.8 


33.47 


16.98 


old 


P 


MCGC3 


00:30:27.30 


+41:36:20.4 


31.88 


16.31 


old 


P 


B514 


00:31:09.90 


+37:53:59.7 


55.39 


15.76 


> 10 Gyr 


G 


MCGC5 


00:35:59.73 


+35:41:03.8 


78.73 


16.09 


old 


P 


B298 


00:38:00.23 


+40:43:55.9 


14.28 


16.59 


old 


C 


H12 


00:38:03.85 


+37:44:00.6 


50.03 


16.47 






B019 


00:40:52.52 


+41:18:53.4 


4.84 


14.93 


old 


c 


B020 


00:40:55.26 


+41:41:25.2 


7.42 


14.91 


interm / old 


P,C 


B023 


00:41:01.18 


+41:13:45.7 


4.46 


14.22 


old 


P,C 


EXT8 


00:53:14.51 


+41:33:24.7 


27.27 


15.54" 






B517 


00:59:59.91 


+41:54:06.6 


45.08 


16.08 







Notes: r p refers to the projected radius from the center of the galaxy. 
" deri ved from ugriz photometry of Peacock et al. (2010) with the transformation equation of Jes ter et al J 
f2005b . 

P: age estimates from lPeacock et al. (2010); C: from Caldwell et al., (20091); G: from lGalleti et all <2005l) . 



We show the spatial distribution of our sample GCs and all the co nfirmed GCs f rom RBC v.4 in 
Figure Q] The large ellipse is the M31 disk/halo boundary as defined bv iRacinel dl99ll) . Note that most 
of our sample are located in the halo of M31 except B019 and B023, which are very close to each other 
with a distance of ~ 5.5 arcmin. That's to say, most of GCs in our sample are halo GCs in M31, which 
can help us to access the nature of galaxy halo with these clusters. 
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Fig. 1 Spatial distribution of our sample GCs (blue filled circles) and all the confirmed GCs 
from RBC v. 4 (open circles). The large ellipse is the M31 disk/halo boundary as defined by 
lRacind(ll99ll) . 



3 OBSERVATIONS AND DATA REDUCTION 

Our Low-resolution spectroscopic observations were all taken at the 2. 16m optical telescope at Xinglong 
Site, which belongs to National Astronomical Observatories, Chinese Academy of Sciences (NAOC), 
from 10th to 13th September 2010. An OMR (Optomechanics Research Inc.) spectrograph and a PI 
1340x400 CCD detector were used during this run with a dispersion of 200 A mm 1 , 4.8 A pixel 1 , 
and a 3.0 "slit. Exposures of 3 x 1800 seconds were taken with seeing typically ~ 2.5 ". Our spectra 
cover the wavelength range of 3500 — 8100 A at 4 A resolution. All our spectra have S/N > 40. 

The spectroscopic data were reduced following standard procedures using the NOAO Image 
Reduction and Analysis Facility (IRAF, version 2.11) software package. The CCD reduction includes 



Spectroscopic Study of M31 Halo GCs 



5 



bias and flat-field correction, as well as cosmic -ray removal. Wavelength calibration was performed 
based on helium/argon lamps exposed at both the beginning and the end of the observations each night. 
Flux calibration of all spectra was perfo rmed based on observations of at least two of the KPNO spec- 
tral standard stars dMassev et all 1 1988b per night. The atmospheric extinction was corrected for using 
the mean extinction coefficients measured for Xinglong by the Beijing-Arizona-Taiwan-Connecticut 
(BATC) multicolor sky survey (H. J. Yan 1995, priv. comm.). 

Before we measure the Lick absorption line indices, the heliocentric radial velocities V r were ob- 
tained by fitting the abosorption lines of our spectra with the templates in various radial velocities. The 
typical internal velocity errors on a single measure is ~ 20 km s _1 . The estimated radial velocities V r 
with the associated uncertainties (Col. 2) are listed in Table [2] The published radial velocities V r (Col. 
3) are also listed for comparisons. It can be seen that our measurements agree well with those listed in 
RBC v.4. At least, we can not see significant differences between our measurements and the published 
values. 

Table 2 The radial velocities V r of our sample GCs as well as 
the previous results. 



ID 


This work 


RBC v.4 


MCGC2 


-586.87 ±39.98 




MCGC3 


-416.46 ± 14.01 




B514 


-429.42 ± 20.24 


-458 ± 23 


MCGC5 


-417.55 ± 25.03 




B298 


-648.50 ± 16.67 


-539 ± 12 


H12 


-412.51 ±33.05 




B019 


-149.83 ± 22.91 


-224 ±2 


B020 


-231.87 ±26.48 


-351 ± 1 


B023 


-348.44 ± 21.30 


-451 ±5 


EXT8 


-104.55 ± 7.32 


-154 ± 30 


B517 


-267.47 ± 20.73 


-272 ± 54 



Similar to iGalleti et al.l d2005l 120061) and ICaldwell et all (I2009I) . we plotted the radial velocity V r 
(corrected for the systemic velocity of M31) versus the projected distance along the major axis (X) in 
Figure [2] The left panel is for all the confirmed clusters while the right panel is for the halo clusters 
which are defined in Figure Q] The small points are the published m easurements from RBC v.4 while 
the filled circles with errors are the measurements in our work. Since ICarignan et all d2006l) calculated 
the HI rotation curve of M3 1 out to ~ 35 kpc with the observations results of Effelsberg and Green Bank 
100 m telescopes, the HI rotation curve of M3 1 galaxy were over plotted in figure[2]with the continuous 
line. It can be seen that both the halo clusters and most of our sample clusters do not follow the disk 
mean velocity curve very well, especially for B514, MCGC5, H12 and B517, suggesting that they do 
not have kinematic association with the star forming young disk of M3 1 . 

Subsequently, all the spectra were shifted to the zero radial velocity and degraded to the wave- 
length dependent Lick reso lution with a variable-width Gaussian kernel following the definition of 
IWorthev & Ottavianil (1 1997b . i.e. 1 1.5 A at 4000 A, 9.2 A at 4400 A, 8.4 A at 4900 A, 8.4 A at 5400 
A, 9.8 A at 6000 A. T hus, we measured all the 25 types of Lick indices stric tly by using the parame- 
ters and formulas from lWorthev et all ( 1 19941) and lWorthev & Ottavianil d 19971). The unce rtainty of each 
index was estimated based on the analytic formulae (111— (18) of ICardiel et alj dl998b . All the Lick 
absorption line indices measurements and ler errors are listed in Table[3] 

As an example, Figure [3] shows the reduced spectroscopy of our sample GC B023, with all the 
Lick absorption line indices bandpasses marked. The spectrum has been degraded and shifted to the 
zero radial velocity as described above. Actually, from the definitions of line indices of CN1 and CN2 
dWorthev et al.l[l994l;IWorthev & Ottavianil[r997l) . we find that the index bandpasses of them are totally 
the same and the only difference is the pseudocontinua coverage. 



6 



Fan et al. 




Fig. 2 The radial velocity V r (corrected for the systemic velocity of M31) as a function 
of the projected distance along the major axis (X) in arcmin for all the confirmed clusters 
(left) and the halo clu sters (right). The solid line is the HI rotation curve of the galaxy from 
ICarignan et alJ d2006l) . The filled circles with errors are the GCs from our sample while the 
small points are the velocity from RBC v.4 catalogue. It is easy to find out that both the 
halo clusters as well as most of our sample clusters are distinct from HI rotation curve of the 
galaxy, implying that they do not have kinematic association with the star forming young disk 
of the galaxy. 



A simple way t o estimate the metall icity is by calculating it from the combination of absorption line 
indices Mg and Fe. Idalleti et al l (l2009h provide the method to measure the metallicity from [MgFe], 
which is defined as [MgFe] = >/Mg6(Fe), with (Fe) = (Fe5270+Fe5335)/2. Thus, the metallicity can 
be derived from the formula below, 

[Fc/H] [Mg Fc] = -2.563 + 1.119[MgFc]-0.106[MgFe] 2 ± 0.15. (1) 

The uncertainty of the [Fc/H][ MgFe ] was estimated with the equation in the following, 

4e/H] = L119 HigFe] + 4 x 0.106 2 [MgFe] 2 af MgFc] . (2) 

All metallicity [Fc/H][ MgF oi derived from [MgFe] and the associated uncertainty determinations are 
listed in Col. (2) of Table |4| It is obvious that all the metallicity derived from the line index [MgFe] 
agree well with those from the model fitting method. 

4 FITTING, ANALYSIS AND RESULTS 
4.1 Model description 

iThomas et al.l d2003l) provided stellar population models including Lick absorption line indices for var- 
ious elemental-abundance ratios, covering ages from 1 to 15 Gyr and metallicities from 1/200 to 3.5 x 
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Indices 


MCGC2 


MCGC3 


B514 


MCGC5 


B298 


H12 


B019 


B020 


B023 


EXT08 


B517 




—2.297 


3.038 


2.526 


3.450 


4.520 


3.331 


1.640 


1.592 





719 


4 


054 


2.797 


error 


0.351 


0.270 


0.315 


0.294 


0.532 


0.249 


0.309 


0.347 





228 





322 


0.258 




0.154 


2.258 


1.916 


2.839 


2.775 


1.977 


0.661 


1.136 





589 


2 


497 


—0.358 


error 


0.237 


0.260 


0.256 


0.307 


0.362 


0.213 


0.233 


0.256 





174 





215 


0.204 


CNl (mag) 


0.037 


-0.084 


-0.095 


-0.137 


-0.104 


-0.069 


-0.039 


-0.056 





031 


-0.142 


-0.097 


error 


0.011 


0.008 


0.007 


0.009 


0.011 


0.007 


0.007 


0.010 





005 





008 


0.006 


CN2 (mag) 


0.097 


0.015 


-0.013 


-0.060 


0.029 


-0.023 


0.033 


0.013 





082 


-0.066 


-0.022 


error 


0.015 


0.022 


0.015 


0.019 


0.024 


0.017 


0.014 


0.012 





011 





020 


0.011 


Ca4227 (A) 


0.358 


0.383 


0.000 


0.119 


0.340 


0.582 


0.388 


0.170 





230 





127 


0.085 


error 


0.219 


0.093 


0.040 


0.110 


0.118 


0.148 


0.097 


0.076 





074 





047 


0.156 


G4300 (A) 


1.416 


1.577 


1.381 


3.088 


0.954 


1.221 


3.117 


2.871 


3 


534 





272 


1.901 


error 


0.482 


0.323 


0.225 


0.396 


0.254 


0.393 


0.319 


0.448 





414 





225 


0.172 


H 7A (A) 


-0.747 


1.064 


0.954 


-1.626 


2.603 


-0.834 


-4.478 


-6.141 


1 


5.271 


1 


746 


-0.215 


error 


0.489 


0.277 


0.249 


0.284 


0.206 


0.339 


0.394 


0.381 





400 





288 


0.239 


H7P (A) 


1.747 


1.364 


1.307 


0.339 


1.346 


1.022 


-0.321 


-0.788 


-0.056 


1 


924 


0.923 


error 


0.254 


0.199 


0.167 


0.252 


0.133 


0.257 


0.157 


0.179 





168 





181 


0.179 


Fe4383 (A) 


-0.202 


—0.202 


-0.520 


—0.654 


0.240 


0.816 


1.957 


2.439 


3 


115 





357 


0.244 


error 


0.543 


0.406 


0.286 


0.240 


0.566 


0.241 


0.335 


0.375 





379 





199 


0.301 


(~!a4455 <h\ 


0.321 


0.399 


0.055 


0.243 


1.309 


0.391 


0.530 


0.291 





601 





118 


0.810 


error 


0.252 


0.224 


0.078 


0.078 


0.260 


0.121 


0.183 


0.099 





127 





030 


0.219 


Fe4531 (A) 


0.169 


0.300 


0.287 


0.823 


0.298 


—0.642 


1.929 


1.711 


1 

1 







281 


2.319 


error 


0.545 


0.212 


0.120 


0.223 


0.208 


0.167 


0.232 


0.348 





155 





111 


0.237 


J_ C±UUO I, jV J 


— 1.582 


—0.848 


2.038 


— 1.400 


0.662 


—2.585 


2.566 


1.162 





848 


-0.498 


1.024 


error 


0.405 


0.289 


0.251 


0.336 


0.433 


0.335 


0.347 


0.181 





214 





080 


0.311 


UB (A) 


2.212 


1.794 


2.250 


2.269 


2.308 


3.082 


2.037 


1.791 


1 


526 


2 


583 


3.148 


error 


0.187 


0.209 


0.211 


0.243 


0.254 


0.305 


0.336 


0.242 





214 





219 


0.290 


Fpsni cA'i 


—0.557 


1.044 


1.506 


1.305 


—0.388 


2.170 


3.209 


3.281 


2 


242 





342 


1.903 


error 


0.386 


0.216 


0.144 


0.198 


0.137 


0.359 


0.384 


0.304 





167 





076 


0.271 


Me:l (mas) 


0.054 


0.004 


-0.007 


-0.003 


0.011 


0.014 


0.034 


0.017 





032 





015 


0.001 


error 


0.005 


0.001 


0.001 


0.002 


0.002 


0.002 


0.002 


0.003 





002 





001 


0.002 


Mg2 (mag) 


0.053 


0.027 


0.028 


0.027 


0.033 


0.036 


0.121 


0.102 





119 





004 


0.044 


error 


0.004 


0.004 


0.003 


0.004 


0.002 


0.002 


0.008 


0.006 





006 





002 


0.004 


Mgb (A) 


0.506 


0.415 


0.518 


1.370 


0.109 


0.731 


2.392 


2.178 


1 


927 





088 


0.962 


error 


0.165 


0.140 


0.098 


0.156 


0.093 


0.105 


0.241 


0.214 





163 





063 


0.164 


Fe5270 (A) 


0.202 


0.409 


0.919 


0.339 


—0.185 


0.148 


1.345 


1.900 


1 


526 





225 


1.618 


error 


0.270 


0.086 


0.140 


0.111 


0.137 


0.103 


0.184 


0.252 





205 





043 


0.194 


Fe5335 (A) 


-0.391 


0.529 


0.166 


1.082 


0.703 


-0.221 


1.014 


0.979 


1 


187 





439 


0.730 


error 


0.217 


0.093 


0.091 


0.187 


0.125 


0.101 


0.190 


0.195 





177 





074 


0.177 


Fe5406 (A) 


0.481 


-0.388 


-0.061 


0.282 


-0.235 


0.234 


0.993 


0.557 





730 





224 


0.108 


error 


0.268 


0.193 


0.091 


0.108 


0.144 


0.115 


0.166 


0.092 





113 





059 


0.101 


Fe5709 (A) 


-0.045 


0.021 


0.008 


0.325 


-0.416 


0.286 


0.351 


0.021 





483 





050 


0.025 


error 


0.082 


0.041 


0.020 


0.061 


0.103 


0.107 


0.075 


0.078 





087 





028 


0.111 


Fe5782 (A) 


0.334 


0.142 


0.177 


0.239 


0.275 


-0.143 


0.311 


0.186 





512 





103 


-0.128 


error 


0.105 


0.028 


0.053 


0.062 


0.094 


0.077 


0.064 


0.053 





097 





029 


0.086 


NaD (A) 


1.025 


1.492 


1.175 


1.446 


1.663 


1.559 


3.491 


2.490 


3 


642 





744 


0.215 


error 


0.113 


0.121 


0.094 


0.164 


0.192 


0.150 


0.399 


0.254 





380 





083 


0.070 


TiOl (mag) 


0.038 


0.004 


0.007 


0.015 


-0.008 


0.009 


0.036 


0.024 





046 





012 


0.002 


error 


0.003 


0.002 


0.002 


0.001 


0.003 


0.002 


0.001 


0.002 





001 





001 


0.004 


Ti02 (mag) 


-0.009 


0.011 


0.002 


0.016 


0.019 


-0.024 


0.062 


0.050 





064 





009 


0.015 


error 


0.003 


0.001 


0.001 


0.001 


0.002 


0.002 


0.001 


0.001 





001 





001 


0.002 



solar abundance. These models are ba s ed on the s tandard models of iMarastonld 19981). wit h input stellar 
evolutionary tr acks from Cassisi et al.l d 1997b and lBono et al. I (119971) and a ISalpeterl(ll955h stellar initial 
mass function jThomas et al. ( 2004 ) improved the models by including higher-order Balmer absorption- 
line indices. They found that these Balmer indices are very sensitive to changes in the a/Fe ratio for su- 
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500 



400 - 




300 — 



200 



4000 



4500 



5000 5500 
Wavelength (A) 



6000 



6500 



Fig. 3 Spectrum of G C B023 in our sample , wit h the index bandpasses of al l the absorption 
Lick indices defined in lWorthev et al.l d 19941) and lWorthev & Ottavianil d 19971) marked. As we 
can see, the index bandpasses for CN1 and CN2 are the same and the only difference is the 
pseudocontinua coverage. 



perso lar metallicities. The lat est stellar popu l ation m ode l for Lick absorption- line indices dThomas et all 
l2010h is an improvement on iThomas et al. I d2003l) and IThomas et all (|2004 . They were derived from 
the MILES stellar library, which provides a higher spectral resolution appropriate for MILES and Sloan 
Digital Sky Survey (SDSS) spectroscopy, as well as flux calibration. The models cover ages from 0.1 
to 15 Gyr, [Fe/H] from —2.25 to 0.67 dex, an d [a/Fe] from —0.3 t o 0.5 dex. In our work, we fitted 
our absorption indices based on the models o f Thomas et al.l d2010l) . by using the two sets of stellar 
evolutionary tracks provided, i.e. JCassisi et al.l (1 1997 ) and Padova. 



4.2 Fitting with stellar population models and the results 

As lCaldwelletaild2009l) demonstrated that the x 2 ~ minimization method for many diagnostic lines are 
more reliable to extract the ages than the two absorption line indices diagram plot method. Furthermore, 
we have measured 25 different types of Lick line indices listed in Table[3j all of w hich were used for the 
fitting procedure, then the results should be much more reliable and accurate. Since lThomas et "all d2010l) 
provide only 20 ages, 6 [Fe/H] values, and 4 [a/Fe], it is necessary to interpolate the original models to 
the higher-resolution models for our needs. We carried out the cubic spline interpolations, using equal 
step lengths, to obtain a grid of 150 ages from 0.1 to 15 Gyr, 31 [Fe/H] values from —2.25 to 0.67 dex, 
and 51 [a/Fe] from —0.3 to 0.5 dex, which makes the model more accurate and more helpful for our 
following statistics. Therefore, the ages (t), metallicities [Fe/H], and [a/Fe] were determined at the same 
time by comparing the interpolated stellar population models with the spectral-energy distributions from 
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Tabl e 4 The y 2 —minimizatio n fitting re s ults using 
iThomas et ail d20ld) models with ICassisi et ail d 19971) and 
Padova stellar tracks, respectively. 



ID [Fe/H] [Mj . F oi [Fe/H] C a S Agc C a S (Gyr) [a/Fe] C a S [Fe/H] Pad Age Pad (Gyr) [a/Fe 



MCGC2 


-2 


32 


± 





28 


-1 


ro+0.28 
JO -0.18 


13 


M+0.40 
ou -0.50 





46« 


02 
7li 


-1 


441°, 


19 

18 


13.60 +0 /oo 





r n +0.00 
du -0.80 


MCGC3 


-2 


09 


± 





13 


-1 


nn+0.09 
8U -0.18 


13 


fi0 +0.50 
OU -2.30 







(10 
47 


-1 


80l° 


IS 
IS 


13.60l°j° 





cn+O.OO 
OU -0.50 


B514 


-2 


00 


± 





12 


-1 


oQ+0.09 
8y -0.18 


13 


fif) +0.40 
DU -0.60 





50±g 


00 

OS 


-1 


891°, 


IS 
IS 


13.60+2TO 





r n +0.00 
du -0.35 


MCGC5 


-1 


56 


± 





18 


-1 


r-o+0.18 
do -0.18 


13 


60 +1.40 
ou -0.60 





501° 


00 
:>,r, 


-1 


441°, 


00 
IS 


13.60liSg 





r n +0.00 
du -0.29 


B298 


-2 


38 


± 





13 


-2 


A7+0.18 
-0.09 


13 


fif) +0.50 





50l° 


00 

Br, 


-2 


071°, 


IX 
18 


13.60l°;|° 





cn+o-oo 

du -0.62 


H12 


-2 


38 


± 





24 


-1 


sn +0.18 


13 


fif) +1.40 





50l° 


00 
20 


-1 


nig 


IS 
18 


14.801°^ 





cn+o.oo 

du -0.20 


B019 


-0 


08 


± 





26 


-0 


7A +0.10 
'^-0.20 


13 


m+0.80 
du -2.70 





481° 


02 
33 


-0 


531°, 


10 
10 


7 7A+6.30 
' ■ ' u -0.80 





,,0+0.02 
48 -0.36 


B020 


-0 


01 


± 





26 


-0 


94+0.20 


13 


yn+0.30 
' U -0.50 





44l° 


00 
41 


-0 


941°, 


10 
10 


13.60l° j° 





cn+OOO 
du -0.35 


B023 


-1 


03 


± 





21 


-0 


0,1+0.10 
84 -0.10 


13 


fi0 +0.50 
ou -0.60 





34l° 


10 
07 


-0 


741°, 


10 
20 


13.201°^ 





4f) +0.10 
^ u -0.34 


EXT8 


-2 


38 


± 





00 


-2 


ny+0.09 
u '-0.18 


13 


m+0.60 
JU -4.30 







00 
SO 

r>() 
00 


-2 


071°, 


Oil 
IS 


13.50li;« 





r n +0.00 
du -0.80 


B517 


-1 


40 


± 





20 


-1 


cq+0.09 
JJ -0.18 


13 


fi0 +1.40 





ootg 


-1 


53±° 


l.s 
09 


13.60l?; 6 7 ° 





m +0.50 
uu -0.30 



Notes: Cas: fitting with the model evolutionary tracks of Cassisi e t all i!997l) ; Pad: fitting with the model 
evolutionary tracks of Padova, respectively. 



our photometry by employing the x 2 — minimization method, i.e., 



xiin(^ Fe / H > a / Fe ) = min 



E 



(3) 



where L™ od (t , Fc/H, a/Fe) is the i Lick line index in the stellar population model for age t, metallicity 
[Fe/H], and [a/Fe], while L^ bs represents the observed Lick line indices from our measurements and 
the errors estimated in our fitting are given as follows, 

a i = O'obs,, + ^mod.i- ( 4 ) 

Here, cr n hs.i; is the ob servational uncertainty and <r mo d,i is the uncertainty associated with the models of 
IThomas et al.l d2010l) . We combined the two uncertainties together in our fitting. 

From TablelU we found that either the [Fe/H] or the ages or the [a/Fe] derived from the lCassisi et al.l 

(Il997h or from Padova tracks of the models are basically the same, suggesting that our fitting results are 
consistent with each other. Moreover, the ages constrained in our work are in good agreement with those 
previous work in Table|4j implying that our fitting method is reliable. Besides, it is worth noting that all 
of our sample GCs, most of which are located in the galaxy halo, are older than 10 Gyr, indicating that 
these clusters formed at the very beginning of the galaxy formation. We also find that the metallicity de- 
rived from the absorption-line index [MgFe] (in C ol. 2) consists with m etallicity fitted with two different 
tracks of the model (in Cols. 3 and 6). Previouslv. lGalleti et al.l d2005l) estimated the metallicity of B514 
with [Fe/H]= -1.8 ± 0.3, and RBC v.4 lists the metallicity for B514 [Fe/H]= -2.06 ± 0.16, for B298 
[Fe/H]= -1.78 ± 0.22, for B019 [Fe/H]= -0.74 ± 0.15, for B020 [Fe/H]= -0.83 ± 0.07, for B023 
[Fe/H]= —0.91 ± 0.14. All these previous measurements are in good agreement with our estimates in 
TablelU indicating our method and r esults are reliable. Fo r consistency, in the su bsequent analysis, we 
adopted the ages and metallicity from lThomas et al] d2010h + ICassisi et al.l d 19971) track of the predictive 
model. 

4.3 Metallicity Properties of Outer Halo 



The metal abundance is one of the most important properties of star clusters to understand the formation 
and enrichment processes of their host galaxy. For instance, the halo stars and clusters should feature 
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large-scale metallicity gradients if the enrichment timescale is shorter than the collapse time, which 
may be due to the galaxy formation as a consequence of a monolithic, dissipative, and rapid coll apse 
of a single mas sive, nearly spherical , spinning gas cloud (lEggen et all ll962tlBarmbv et"al l l2000l) . On 
the other hand. lSearle & Zinnld 19781) presented a chaotic scheme for early galactic evolution, when the 
loosely bound pre-enriched fragments merge with the protogalaxy during a very long period of time, in 
which case a more homogeneous metallicity distribution should develop. Most galaxies are thought to 
have formed through a combination of these scenarios. 

A lot of previous work attempted to find the clues of formation and evol ution for M31 galax y 
through studying the metallicity distr i bution of its globu l ar clu st er system . jHuchra et al.l dl99ll) : 
lAshman & Birdl dl993l) : IBarmbv et all d2000t) ; IPerrett et all d2002l) ; iFan et al.l (12008b found that the 
metal-rich GCs are statistically more concentrated toward the center of the galaxy, while their metal- 
poor counterparts are more spatially extended in the halo. Furthermore, there are also many work to fin d 
out whether or not a radial metallicity gradient exists for M3 1 star cluster system. Ivan denBerghld!969l) : 
iHuchraet al.l d 19821) showed that there is little or no evidence for a general radial metallicity gradient fo r 
GCs with i n a rad ius of 50 arcmin. However, studies including lHuchra et al.l (1 1 99 1 b : IPerrett et al.l d2002l) ; 
IFan et all d2008l) support the possible existenc e of a radial metallic ity gradient for the metal-poor M3 1 
GCs, although the slope is not very significant. IPerrett et al.ld2002l) suggest that the gradients is —0.017 
and — 0.015 dex arcmin -1 for the full sample and inner metal-poor clusters. More recently. IFan et al.l 
d2008l) found that the slope is —0.006 and —0.007 dex arcmin -1 for the metal-poor subsample and 
whole sample while the slope approches zero for the metal-rich subsample. Nevertheless, all these stud- 
ies are based on GCs that are located relatively close to the center of the galaxy, usually at projected 
radii of less than 100 arcmin. In our work, we extended the radial coverage to a radius of r p ~ 117 
kpc, which corresponds to ~ 510 arcmin, to check if the previous findings are correct at a much larger 
distance from the galaxy center. 

For the purpose of better investigating the metallicity distribution/spatial gradient, we enlarged the 
metallicity samp le by merging the metallicity of ou r meas u rements with the published spectroscopi c 
metallicity fromlHuchra et al.1 dl99ll): IBarmbv et all d2000l) ; IPerrett et alj d2002l): iGalleti et al.l d2009l) ; 
ICaldwell et all d201ll) together with those from CMD fittings iMackev et alj d2006l 120071 12010|) . For 
the published data, if the metallicity from different work overlapped with the other, the small smaller 
associated uncertainty data will superseded the larger one and the spectroscopic data will superseded 
the one derived from CMD fitting. In total, we have a metallicity sample of 384 entries. 

Figure [4] shows the metallicity as a function of projected radius from the galaxy center for all 
M3 1 confirmed clusters (Top) and the halo clusters only (Bottom) in the unit of kpc. In the top panel, 
o pen triangles with err o rbars represent the sp e ctroscopic metallicit i es from the published measurements 
of iHuchra etail dl99ll) : lBarmbv et alj d2000l) : IPerrett et al.l d2002l);lGalleti et alj d2009l) : ICaldwell et al.l 



d201 1 ) as well as the metallicities from CMD fittings IMackev et al. I (120061 12007L 120 1 Oh while the green 



filled triangles with errorbars are our spectroscopic measurements. The solid line is a linear fit to all 
the data points, with a slope of —0.028 ± 0.001 dex kpc -1 , resp onding to —0.007 dex arcmin" 1 . Th e 
fit results are similar to those given by the previous works (see, IPerrett et all 120021; IFan et al L l2008h . 
which are based on the cluster sample within a projected radii ?' p < 100 arcmin (~ 23 kpc). Thus, our 
work tentatively supports the notion that a radial metallicity gradient may exist out to a projected radius 
of ~ 117 kpc by merging the published metallicities. In other words, we updated the results with the 
new sample extended to M31's most remote outer halo. Since the aim of our work is to study the nature 
of M31 halo, we would like to foucs on the metallicity gradient of halo clusters. In the bottom panel, 
we only plot the halo clusters, which are defined in Figure Q] A least-squares fitting yields the slope of 
—0.018 ± 0.001 dex kpc -1 . Therefore, it can be seen that the metallicity gradient seems to exist for the 
halo clusters, although it is not significant. 

Furthermore, It is noted that in Figure|4j the metallicity gradient for the clusters located in the outer 
halo with r p > 25 kpc is not significant. Thus, we show the metallicity as a function of projected radius 
for only the outer halo clusters with r p > 25 kpc in Figure|5] A least-squares linear fitting for all the data 
shows the slope is —0.010 ± 0.002 dex kpc -1 (the solid line). However, if G001 is excluded, the slope 
turns out to be —0.004 ± 0.002 dex kpc -1 (the red dashed line), which is much shallower than that in 





±, 20 40 60 80 100 120 



r p (kpc) 



Fig. 4 [Fe/H] versus projected radius from the galaxy center for M31 GCs. The solid line 
refers to a linear fit to all the data. Top: All the confirmed clusters. The open triangles with 
errorbars represent published metallicities while the green filled triangles with errorbars are 
our measurements. Bottom: Halo clusters only. All the clusters are marked with filled circles. 



Figure|4] It m ay suggest that the m etallicity gradient is not significant for the outer halo clusters in M3 1 . 
Very recently. iHuxor et al.l d201 lb investigated the metallicity grad i ent for 15 halo CG s to r p =117 kpc 
with the metallicity derived from the CMD fittings iMackev etail d2006l 120071 12010|) and the authors 
found that the metallicity gradient becomes not significant if one halo G C H14 is excluded in their 
Figure 6. We found that our result is consistent with the previous finding o f lHuxoret al.ld201ll) . 

5 DISCUSSIONS AND SUMMARY 

In our work, we carried out the spectroscopic observations of 1 1 confirmed globular clusters of M3 1 
with the OMR spectrograph and a PI 1340x400 CCD detector on 2.16 m telescope at Xinglong site of 
NAOC from 10th to 13th September of 2010. Since our aim is to study the nature of the halo of M31, 
we selected the bright confirmed clusters, 9 of which are located in the halo, out to a projected radius of 
78.75 kpc from the galactic center. 

For all ou r sample clust e rs, we measured all types of the Lick absorption-line indices (see the 
definitions in, IWorthev et al.L 1 19941 IWorthev & Ottavianii 1 19971) as well as the radial velocities. We 
found that most GCs of our sample are distinct from the HI rotation curve of M3 1 galaxy, especially 
for B514, MCGC5, H12 and B517, suggesting that most of our sample clusters do not have kinematic 
association with the star forming young disk of the galaxy. 
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> -2 



-3 



fit to all data 
excluding G001 
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20 



AO 



60 



30 



00 



20 
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Fig. 5 [Fe/H] from versus projected radius from the galaxy center for the outer halo GCs, 
whose r p > 25 kpc from the center of the galaxy. The slope of the fitting for all the data is 



-0.010 ± 0.002 dex kpc" 1 . 
0.002 dex kpc- 1 . 



However, if G001 is excluded, the slope turns out to be —0.004 ± 



Since ICaldwell et ail (120091) demonstrated that the % 2 — minimization method for many diagnostic 
lines are more reliable for extracting the ages than the line indices diagram plot, in our wok we ap- 
plied the \ 2 — minim ization method to fit the line indices with the updated stellar population model 
iThomas et al.l (12010b with two different tracks Cassisi and Padova, separately. The fitting results show 
that all our sample clusters are older than 10 Gyr and most of them are metal-poor (—0.91 < [Fe/H] 
< -2.38 dex). 

In order to enlarge our sample, we merged the spectroscopic metallicity of our work with the pre- 
viously published ones, extending the cluster sample out to a projected radius of 117 kpc from the 
galaxy's center. We found the metallicity gradient for all the confirmed clusters exists with a slope of 
-0.028 ± 0.001 dex kpc -1 . However, the slope turns to be -0.018 ± 0.001 dex kpc" 1 for all the halo 
clusters, which is much shallower. If we only consider the outer halo clusters with r p > 25 kpc, the 
slope becomes —0.010 ± 0.002 dex kpc -1 and if one cluster G001 is excluded from the outer halo sam- 
ple, the slope even could be —0.004 ± 0.002 dex kpc -1 . Thus we conclude that metallicity gradient for 
M31 outer halo clusters is not significant, which agrees well with the previous findings. This result may 
imply that for M31 galaxy formation, the "rapid collapsing" scenario is supported inside the inner halo 
while the "fragments merging" scenario is proved in the outer halo of the galaxy beyond 25 kpc from 
the center. It seems that the combination of the two scenarios could best explain the galaxy formation. 
However, we still need more observations and further study to figure it out. 

Besides, it is interesting to note that the halo of M3 1 galaxy might be divided into two parts (by 
combining the Huxor data): inner halo and outer halo from our study. The nature of the two parts of 
halo seem to be different in terms of metallicity gradients of the star clusters, which may be due to the 
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different formation mechanisms of the two parts. Just like the Milky Way halo from the SDSS/SEGUE 
data, the Milky Way halo could be divided into two parts with different metallicity properties based on 
the observations of a large sample halo stars. Therefore, it seems that M31 galaxy and our Galaxy have 
more similarities than we expected. However, more observational data is required for the further study 
in the future. 
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